Analytical method of determining the solution of an object with edges under applied load

ABSTRACT

An analytical method of determining the solution of an object with a plurality of planar surfaces under applied load includes the steps of assigning the load applied on each of the plurality of planar surfaces in the three dimensional directions respectively, determining the equivalent load applied on each of the plurality of planar surfaces in the three dimensional directions from the applied loads with the stress boundary conditions, and determining the stress and strain fields from the equivalent load applied on each of the plurality of planar surfaces.

TECHNICAL FIELD

The present invention relates to an analytical method of determining the solution (deformation and stress) of an object, and particularly, although not exclusively, to an analytical method of determining the deformation and stress of an object with a plurality of planar surfaces under applied load.

BACKGROUND

Contact mechanics is of particular importance in determining the physical behaviour of a mechanical component such as rail-wheel systems, roller bearings, and cutting tools, where the contact or loaded areas of these components are not far from the edge surface of the contact bodies.

In such situations, the edge surface may affect the internal stress and strain fields, and thus it is impossible to obtain accurate analytical solutions using classical contact mechanic formulae which are derived based on the assumption that the loaded object is an elastic half-space.

SUMMARY OF THE INVENTION

In accordance with a first aspect of the present invention, there is provided an analytical method of determining the solution of an object with a plurality of planar surfaces under applied load, comprising the steps of:

A. assigning the load applied on each of the plurality of planar surfaces in the three dimensional directions respectively; and

B. determining the equivalent load applied on each of the plurality of planar surfaces in the three dimensional directions from the applied load with the stress boundary conditions.

C. determining the stress and strain fields from the equivalent load applied on each of the plurality of planar surfaces.

In an embodiment of the first aspect, each of the plurality of planar surfaces is divided into a plurality of uniformly loaded sections each having a centre respectively.

In an embodiment of the first aspect, the applied load includes the accumulation of uniform loads applied on each uniformly loaded section.

In an embodiment of the first aspect, the equivalent load includes the accumulation of uniform equivalent loads applied on each uniformly loaded section.

In an embodiment of the first aspect, the induced stress includes the accumulation of stresses induced on the centre of each uniformly loaded section.

In an embodiment of the first aspect, each centre of the uniformly loaded sections on the plurality of planar surfaces X₁-Y₁ and X₂-Y₂ is represented by x_(1i), y_(1i) and x_(2j), y_(2j) respectively.

In an embodiment of the first aspect, the parameter of the applied load applied on each of the plurality of planar surfaces X₁-Y₁ and X₂-Y₂ are represented by P₁, Q_(X1), Q_(Y1) and P₂, Q_(X2), Q_(Y2) respectively.

In an embodiment of the first aspect, P₁, Q_(X1), Q_(Y1) and P₂, Q_(X2), Q_(Y2) are each represented in the form of the accumulation of uniform loads applied on each uniformly loaded section by the following expression respectively:

${P_{1} = \begin{bmatrix} \left( p_{1} \right)_{1} \\ \left( p_{1} \right)_{2} \\ \vdots \\ \left( p_{1} \right)_{k} \end{bmatrix}};{Q_{X\; 1} = \begin{bmatrix} \left( q_{x\; 1} \right)_{1} \\ \left( q_{x\; 1} \right)_{2} \\ \vdots \\ \left( q_{x\; 1} \right)_{k} \end{bmatrix}};{Q_{Y\; 1} = \begin{bmatrix} \left( q_{y\; 1} \right)_{1} \\ \left( q_{y\; 1} \right)_{2} \\ \vdots \\ \left( q_{y\; 1} \right)_{k} \end{bmatrix}}$ ${P_{2} = \begin{bmatrix} \left( p_{2} \right)_{1} \\ \left( p_{2} \right)_{2} \\ \vdots \\ \left( p_{2} \right)_{l} \end{bmatrix}};{Q_{X\; 2} = \begin{bmatrix} \left( q_{x\; 2} \right)_{1} \\ \left( q_{x\; 2} \right)_{2} \\ \vdots \\ \left( q_{x\; 2} \right)_{l} \end{bmatrix}};{Q_{Y\; 2} = \begin{bmatrix} \left( q_{y\; 2} \right)_{1} \\ \left( q_{y\; 2} \right)_{2} \\ \vdots \\ \left( q_{y\; 2} \right)_{l} \end{bmatrix}}$ where (p₁)_(i)=p₁(x_(1i), y_(1i)), (q_(x1))_(i)=q_(x1)(x_(1i), y_(1i)), (q_(y1))_(i)=q_(y1)(x_(1i), y_(1i)); and (p₂)_(j)=p₂(x_(2j), y_(2j)), (q_(x2))_(j)=q_(x2)(x_(2j), y_(2j)), (q_(y2))_(j)=q_(y2)(x_(2j), y_(2j)).

In an embodiment of the first aspect, the parameter of the equivalent load applied on the plurality of planar surfaces X₁-Y₁ and X₂-Y₂ are each represented by P ₁, Q _(X1), Q _(Y1) and P ₂, Q _(X2), Q _(Y2) respectively.

In an embodiment of the first aspect, P ₁, Q _(X1), Q _(Y1) and P ₂, Q _(X2), Q _(Y2) are each represented in the form of the accumulation of uniform equivalent loads applied on each uniformly loaded section by the following expression respectively:

${{\overset{\_}{P}}_{1} = \begin{bmatrix} \left( {\overset{\_}{p}}_{1} \right)_{1} \\ \left( {\overset{\_}{p}}_{1} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{p}}_{1} \right)_{k} \end{bmatrix}};{{\overset{\_}{Q}}_{X\; 1} = \begin{bmatrix} \left( {\overset{\_}{q}}_{x\; 1} \right)_{1} \\ \left( {\overset{\_}{q}}_{x\; 1} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{q}}_{x\; 1} \right)_{k} \end{bmatrix}};{{\overset{\_}{Q}}_{Y\; 1} = \begin{bmatrix} \left( {\overset{\_}{q}}_{y\; 1} \right)_{1} \\ \left( {\overset{\_}{q}}_{y\; 1} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{q}}_{y\; 1} \right)_{k} \end{bmatrix}}$ ${{\overset{\_}{P}}_{2} = \begin{bmatrix} \left( {\overset{\_}{p}}_{2} \right)_{1} \\ \left( {\overset{\_}{p}}_{2} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{p}}_{2} \right)_{l} \end{bmatrix}};{{\overset{\_}{Q}}_{X\; 2} = \begin{bmatrix} \left( {\overset{\_}{q}}_{x\; 2} \right)_{1} \\ \left( {\overset{\_}{q}}_{x\; 2} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{q}}_{x\; 2} \right)_{l} \end{bmatrix}};{{\overset{\_}{Q}}_{Y\; 2} = \begin{bmatrix} \left( {\overset{\_}{q}}_{y\; 2} \right)_{1} \\ \left( {\overset{\_}{q}}_{y\; 2} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{q}}_{y\; 2} \right)_{l} \end{bmatrix}}$ where (p ₁)_(i)=p ₁(x_(1i), y_(1i)), (q _(x1))_(i)=q _(x1)(x_(1i), y_(1i)), (q _(y1))_(i)=q _(y1)(x_(1i), y_(1i)); and (p ₂)_(j)=p ₂(x_(2j), y_(2j)), (q _(x2))_(j)=q _(x2)(x_(2j), y_(2j)), (q _(y2))_(j)=q _(y2)(x_(2j), y_(2j)).

In an embodiment of the first aspect, the parameter of the induced stress induced by the equivalent load on the plurality of planar surfaces X₂-Y₂ and X₁-Y₁ are each represented by S_(Z), S_(X), S_(Y) and T_(X), T_(X), T_(Y) respectively.

In an embodiment of the first aspect, S_(Z), S_(X), S_(Y) and T_(Z), T_(X), T_(Y) are each represented in the form of the accumulation of stresses induced on the centre of each uniformly loaded section by the following expression respectively:

${S_{Z} = \begin{bmatrix} \left( \sigma_{{zz}\; 2} \right)_{1} \\ \left( \sigma_{{zz}\; 2} \right)_{2} \\ \vdots \\ \left( \sigma_{{zz}\; 2} \right)_{l} \end{bmatrix}};{S_{X} = \begin{bmatrix} \left( \tau_{{zx}\; 2} \right)_{1} \\ \left( \tau_{{zx}\; 2} \right)_{2} \\ \vdots \\ \left( \tau_{{zx}\; 2} \right)_{l} \end{bmatrix}};{S_{Y} = \begin{bmatrix} \left( \tau_{{zy}\; 2} \right)_{1} \\ \left( \tau_{{zy}\; 2} \right)_{2} \\ \vdots \\ \left( \tau_{{zy}\; 2} \right)_{l} \end{bmatrix}}$ ${T_{Z} = \begin{bmatrix} \left( \sigma_{{zz}\; 1} \right)_{1} \\ \left( \sigma_{{zz}\; 1} \right)_{2} \\ \vdots \\ \left( \sigma_{{zz}\; 1} \right)_{k} \end{bmatrix}};{T_{X} = \begin{bmatrix} \left( \tau_{{zx}\; 1} \right)_{1} \\ \left( \tau_{{zx}\; 1} \right)_{2} \\ \vdots \\ \left( \tau_{{zx}\; 1} \right)_{k} \end{bmatrix}};{T_{Y} = \begin{bmatrix} \left( \tau_{{zy}\; 1} \right)_{1} \\ \left( \tau_{{zy}\; 1} \right)_{2} \\ \vdots \\ \left( \tau_{{zy}\; 1} \right)_{k} \end{bmatrix}}$ where (σ_(zz2))_(j)=σ_(zz2)(x_(2j), y_(2j)), (τ_(zx2))_(j)=τ_(zx2)(x_(2j), y_(2j)), (τ_(zy2))_(j)=τ_(zy2)(x_(2j), y_(2j)); and (σ_(zz1))_(i)=σ_(zz1)(x_(1i), y_(1i)), (τ_(zx1))_(i)=τ_(zx1)(x_(1i), y_(1i)), (τ_(xy1))_(i)=τ_(zy1)(x_(1i), y_(1i)).

In an embodiment of the first aspect, the stresses (σ_(zz2))_(j), (τ_(zx2))_(j), (τ_(zy2))_(j) on the centre of the jth uniformly loaded section on planar surface X₂-Y₂ by equivalent loads (p _(i))_(i), (q _(x1))_(i), (q _(y1))_(i) on the centre of the uniformly loaded section ith on planar surface X₁-Y₁ are determined by the following expression: σ_(zz2)(x _(2j) ,y _(2j))=m _(zzij) p ₁(x _(1i) ,y _(1i))+m _(zxij) q _(x1)(x _(1i) ,y _(1i))+m _(zyij) q _(y1)(x _(1i) ,y _(1i)) τ_(xz2)(x _(2j) ,y _(2j))=m _(xzij) p ₁(x _(1i) ,y _(1i))+m _(xxij) q _(x1)(x _(1i) ,y _(1i) +m _(xyij) q _(y1)(x _(1i) ,y _(1i)) τ_(zy2)(x _(2j) ,y _(2j))=m _(yzij) p ₁(x _(1i) ,y _(1i))+m _(yxij) q _(x1)(x _(1i) ,y _(1i))+m _(yyij) q _(y1)(x _(1i) ,y _(1i))

In an embodiment of the first aspect, the stresses (σ_(zz1))_(i), (σ_(zx1))_(i), (σ_(zy1))_(i) on the centre of the ith uniformly loaded section on planar surface X₁-Y₁ by equivalent loads (p ₂)_(j), (q _(x2))_(j), (q _(y2))_(j) on the centre of the jth uniformly loaded section on planar surface X₂-Y₂ are determined by the following expression: σ_(zz1)(x _(1i) ,y _(1i))=n _(zxji) p ₂(x _(2j) ,y _(2j))+n _(zxji) q _(x2)(x _(2j) ,y _(2j))+n _(zyji) q _(y2)(x _(2j) ,y _(2j)) τ_(zx1)(x _(1i) ,y _(1i))=n _(xzji) p ₂(x _(2j) ,y _(2j))+n _(xxji) q _(x2)(x _(2j) ,y _(2j))+n _(xyji) q _(y2)(x _(2j) ,y _(2j)) τ_(zy1)(x _(1i) ,y _(1i))=n _(yzji) p ₂(x _(2j) ,y _(2j))+n _(yxji) q _(x2)(x _(2i) ,y _(2j))+n _(yyji) q _(y2)(x _(2j) ,y _(2j))

In an embodiment of the first aspect, the coefficients M_(zz), m_(xz), m_(yz) and n_(zz), n_(xz), n_(yz) are determined based on the classical half-space formulae.

In an embodiment of the first aspect, the stress and strain fields induced by the applied load on the plurality of planar surfaces is determined based on the corresponding equivalent load using the classical half-space formulae.

In accordance with a second aspect of the present invention, there is provided an analytical method of determining the solution of an object with a plurality of planar surfaces under applied load, comprising the step of determining a transformation matrix and in turn providing the equivalent load from the original applied-load.

In an embodiment of the second aspect, the transformation matrix is associated with the physical characteristic of a said object.

In an embodiment of the second aspect, the physical characteristic includes the elastic property of a said object and the relative orientations between the planar surfaces.

In an embodiment of the second aspect, the transformation matrix is associated with the size of grids or meshes of the uniformly loaded section.

In an embodiment of the second aspect, the transformation matrix includes the accumulation of coefficients used in the classical formulae for half-space stress and strain calculations.

In an embodiment of the second aspect, the transformation matrix is represented by a plurality of two dimensional coefficients M_(ZZ), M_(ZX), M_(ZY), M_(XZ), M_(XX), M_(XY), M_(YZ), M_(YX), M_(YY) and N_(ZZ), N_(ZX), N_(ZY), N_(XZ), N_(ZZ), N_(XY), N_(YZ), N_(YX), N_(YY) respectively.

In an embodiment of the second aspect, M_(ZZ), M_(ZX), M_(ZY), M_(XZ), M_(XX), M_(XY), M_(YZ), M_(YX), M_(YY) and N_(ZZ), N_(ZX), N_(ZY), N_(XZ), N_(XX), N_(XY), N_(YZ), N_(YX), N_(YY) are each represented in the form of the accumulation of coefficients of each uniformly loaded section by the following expression respectively:

$M_{ZZ} = \begin{bmatrix} m_{{zz}\; 11} & m_{{zz}\; 21} & \ldots & m_{{zzk}\; 1} \\ m_{{zz}\; 12} & m_{{zz}\; 22} & \ldots & m_{{zz}\; k\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ m_{{zz}\; 1l} & m_{{zz}\; 2l} & \ldots & m_{{zz}\;{kl}} \end{bmatrix}$ $M_{ZX} = \begin{bmatrix} m_{{zx}\; 11} & m_{{zx}\; 21} & \ldots & m_{{zxk}\; 1} \\ m_{{zx}\; 12} & m_{{zx}\; 22} & \ldots & m_{{zx}\; k\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ m_{{zx}\; 1l} & m_{{zx}\; 2l} & \ldots & m_{{zx}\;{kl}} \end{bmatrix}$              ⋮ $M_{YY} = \begin{bmatrix} m_{{yy}\; 11} & m_{{yy}\; 21} & \ldots & m_{{yy}\; k\; 1} \\ m_{{yy}\; 12} & m_{{yy}\; 22} & \ldots & m_{{yy}\; k\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ m_{{yy}\; 1l} & m_{{yy}\; 2l} & \ldots & m_{{yy}\;{kl}} \end{bmatrix}$ $N_{ZZ} = \begin{bmatrix} n_{{zz}\; 11} & n_{{zz}\; 21} & \ldots & n_{{zzl}\; 1} \\ n_{{zz}\; 12} & n_{{zz}\; 22} & \ldots & n_{{zz}\; l\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ n_{{zz}\; 1k} & n_{{zz}\; 2k} & \ldots & n_{{zz}\;{lk}} \end{bmatrix}$ $N_{ZX} = \begin{bmatrix} n_{{zx}\; 11} & n_{{zx}\; 21} & \ldots & n_{{zxl}\; 1} \\ n_{{zx}\; 12} & n_{{zx}\; 22} & \ldots & n_{{zx}\; l\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ n_{{zx}\; 1k} & n_{{zx}\; 2k} & \ldots & n_{{zx}\;{lk}} \end{bmatrix}$              ⋮ $N_{YY} = \begin{bmatrix} n_{{yy}\; 11} & n_{{yy}\; 21} & \ldots & n_{{yy}\; l\; 1} \\ n_{{yy}\; 12} & n_{{yy}\; 22} & \ldots & n_{{yy}\; l\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ n_{{yy}\; 1k} & n_{{yy}\; 2k} & \ldots & n_{{yy}\;{lk}} \end{bmatrix}$

In an embodiment of the second aspect, following step B of the first aspect, P ₁, Q _(X1), Q _(Y1) and P ₂, Q _(X2), Q _(Y2) are determined from the following expression:

$\begin{bmatrix} {\overset{\_}{P}}_{1} \\ {\overset{\_}{Q}}_{X\; 1} \\ {\overset{\_}{Q}}_{Y\; 1} \\ {\overset{\_}{P}}_{2} \\ {\overset{\_}{Q}}_{X\; 2} \\ {\overset{\_}{Q}}_{Y\; 2} \end{bmatrix} = {\begin{bmatrix} I & 0 & 0 & {- N_{ZZ}} & {- N_{ZX}} & {- N_{ZY}} \\ 0 & I & 0 & {- N_{XZ}} & {- N_{XX}} & {- N_{XY}} \\ 0 & 0 & I & {- N_{YZ}} & {- N_{YX}} & {- N_{YY}} \\ {- M_{ZZ}} & {- M_{ZX}} & {- M_{ZY}} & I & 0 & 0 \\ {- M_{XZ}} & {- M_{XX}} & {- M_{XY}} & 0 & I & 0 \\ {- M_{YZ}} & {- M_{YX}} & {- M_{YY}} & 0 & 0 & I \end{bmatrix}^{- 1}\begin{bmatrix} P_{1} \\ Q_{X\; 1} \\ Q_{Y\; 1} \\ P_{2} \\ Q_{X\; 2} \\ Q_{Y\; 2} \end{bmatrix}}$ where the inverse matrix is the transformation matrix.

In an embodiment of the second aspect, the transformation matrix is represented with the following expression:

$\begin{bmatrix} I & 0 & 0 & {- N_{ZZ}} & {- N_{ZX}} & {- N_{ZY}} \\ 0 & I & 0 & {- N_{XZ}} & {- N_{XX}} & {- N_{XY}} \\ 0 & 0 & I & {- N_{YZ}} & {- N_{YX}} & {- N_{YY}} \\ {- M_{ZZ}} & {- M_{ZX}} & {- M_{ZY}} & I & 0 & 0 \\ {- M_{XZ}} & {- M_{XX}} & {- M_{XY}} & 0 & I & 0 \\ {- M_{YZ}} & {- M_{YX}} & {- M_{YY}} & 0 & 0 & I \end{bmatrix}^{- 1}$

BRIEF DESCRIPTION OF THE DRAWINGS

Embodiments of the present invention will now be described, by way of example, with reference to the accompanying drawings in which:

FIG. 1 is an illustration showing a three-dimensional wedge problem;

FIG. 2 is a flowchart showing an example embodiment of the method in the present invention;

FIG. 3a is an illustration showing a wedge with normal and shear loads (p and q);

FIG. 3b is an illustration showing the superposition of two half-spaces with equivalent loads (p, and q) and corresponding induced stresses on the two planar surfaces of the wedge;

FIG. 4a is an illustration showing the X₁-Y₁ plane of the wedge with a plurality of grids;

FIG. 4b is an illustration showing an individual rectangle of the grids;

FIG. 5 is an illustration showing a half-semispherical normal load distribution applied on the top surface of the wedge structure;

FIG. 6 is an illustration showing the applied load on the top surface of the wedge structure in three dimensional directions;

FIG. 7a is an illustration showing the calculated equivalent normal load on X₁-Y₁ plane of a 170° wedge;

FIG. 7b is an illustration showing the calculated equivalent shear load in X₁ direction on X₁-Y₁ plane of a 170° wedge;

FIG. 7c is an illustration showing the calculated equivalent shear load in Y₁ direction on X₁-Y₁ plane of a 170° wedge;

FIG. 7d is an illustration showing the calculated equivalent normal load on X₂-Y₂ plane of a 170° wedge;

FIG. 7e is an illustration showing the calculated equivalent shear load in X₂ direction on X₂-Y₂ plane of a 170° wedge;

FIG. 7f is an illustration showing the calculated equivalent shear load in Y₂ direction on X₂-Y₂ plane of a 170° wedge;

FIG. 8a is an illustration showing the normal stress distribution σ_(xx) in the Y=0 plane of a 170° wedge;

FIG. 8b is an illustration showing the normal stress distribution σ_(yy) in the Y=0 plane of a 170° wedge;

FIG. 8c is an illustration showing the normal stress distribution σ_(zz) in the Y=0 plane of a 170° wedge;

FIG. 8d is an illustration showing the von Mises stress in the Y=0 plane of a 170° wedge;

FIG. 9a is an illustration showing the normal stress distribution σ_(xx) in the Y=0 plane of a 90° wedge;

FIG. 9b is an illustration showing the normal stress distribution σ_(yy) in the Y=0 plane of a 90° wedge;

FIG. 9c is an illustration showing the normal stress distribution σ_(zz) in the Y=0 plane of a 90° wedge;

FIG. 9d is an illustration showing the von Mises stress in the Y=0 plane of a 90° wedge;

FIG. 10a is an illustration showing the normal stress distribution σ_(xx) in the Y=0 plane of a 60° wedge;

FIG. 10b is an illustration showing the normal stress distribution σ_(yy) in the Y=0 plane of a 60° wedge;

FIG. 10c is an illustration showing the normal stress distribution σ_(zz) in the Y=0 plane of a 60° wedge;

FIG. 10d is an illustration showing the von Mises stress in the Y=0 plane of a 60° wedge;

FIG. 11a is an illustration showing the normal stress distribution σ_(xx) in the Y=0 plane of a 135° wedge;

FIG. 11b is an illustration showing the normal stress distribution σ_(yy) in the Y=0 plane of a 135° wedge;

FIG. 11c is an illustration showing the normal stress distribution σ_(zz) in the Y=0 plane of a 135° wedge;

FIG. 12a is an illustration showing the normal stress distribution σ_(xx) in the Y=0 plane of a 75° wedge;

FIG. 12b is an illustration showing the normal stress distribution σ_(yy) in the Y=0 plane of a 75° wedge;

FIG. 12c is an illustration showing the normal stress distribution σ_(zz) in the Y=0 plane of a 75° wedge;

FIG. 13a is an illustration showing the normal stress distribution σ_(zz) in the Y=0 plane of a 90° wedge under two half-semispherical distributed loads on X₁-Y₁ plane and X₂-Y₂ plane;

FIG. 13b is an illustration showing the normal stress distribution σ_(yy) in the Y=0 plane of a 90° wedge under two half-semispherical distributed loads on X₁-Y₁ plane and X₂-Y₂ plane; and

FIG. 13c is an illustration showing the von Mises stress in the Y=0 plane of a 90° wedge under two half-semispherical distributed loads on X₁-Y₁ plane and X₂-Y₂ plane.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

The inventors have, through their own research and experiences, devised that enabling the determination of the induced stress and deformation of an object under load is utterly important in engineering. A semi-infinite (termed as half-space in contact mechanics) body can be solved analytically with classical formulae.

However, 3D objects are generally composed of various planar surfaces or cross sections and how they respond to external loads cannot be directly calculated. The use of alternative methods such as the finite element method (FEM) is thus very common in practice. However, high precision FEM results can only be achieved using very fine meshes, which mean very long computational time and exceedingly large memories.

The inventors have developed an algorithm to derive the stress and strain solutions of three-dimensional (3D) objects. Its concept can be illustrated using a simple 3D object, wedge, which is constructed with two planar surfaces. A wedge under applied loads on its two planes is made equivalent to the superposition of two separated half-spaces under equivalent loads. The equivalent loads are directly related to the original applied-loads of the wedge and can be deduced based on the original stress boundary condition of the wedge (being the original loads on its two planes). Knowing the equivalent loads, the stress and strain of the wedge 10 can be determined by using the classical formulae of half-spaces.

Referring initially to FIG. 1, there is shown geometry of a wedge 10 having a loaded top surface 20, an edge surface 30, and a wedge angle θ of <180°. The loading on the top surface 20 is subjected to a concentrated or any other type of distributed load, and thus give rise to a 3D wedge problem.

The present invention adopts an analytical approach to tackle such 3D elastic wedge problem through the stress analysis of the wedge 10. Both normal and shear stresses on the two planes of the wedge 10 are considered in the analysis. In this approach, such stress analysis is conducted by using matrix operations, and thus explicit solutions may be obtained. Accordingly, the present invention provides a feasible solution to the stress problem of 3D wedge with given loads, i.e. a fixed stress boundary condition applied.

With reference to FIG. 2, there is provided an example embodiment of an analytical method 100 of determining the solution of an object 10 with a plurality of planar surfaces 20, 30 under applied load, comprising the steps of:

-   -   A. assigning given values to the load (p and q) applied on each         of the plurality of planar surfaces 20, 30 in the three         dimensional directions respectively (step 102);     -   B. determining the transformation matrix comprising coefficients         associated with the physical characteristic of a said object and         thus the equivalent load (p and q) applied on each of the         plurality of planar surfaces 20, 30 in the three dimensional         directions from the applied load (p and q) at stress boundary         condition (step 104);     -   C. determining the internal stress field (σ and τ) induced by         the applied load (p and q) applied on each of the plurality of         planar surfaces 20, 30 based on the equivalent load (p and q)         using half-space formulae (step 106);     -   D. determining the internal strain field induced by the applied         load (p and q) applied on each of the plurality of planar         surfaces 20, 30 based on the equivalent load (p and q) using         half-space formulae (step 108).

The inventors of the present invention have derived that the equivalent load p and q of half-spaces can be determined from the applied load p and q and the transformation matrix which is associated with the physical characteristic of the object 10, such as the elastic property of the object 10 and the relative orientations between the planar surfaces 20, 30 (step 104).

Accordingly, once the equivalent normal p and shear q load are determined from the method of the present invention, the stress and strain fields induced in the body 10 can be obtained by using the classical formulae derived based on half-space model (steps 106, 108). Finally, the initiation of mechanical failure taking the form of plastic yielding and deformation of the object 10 due to the applied load p and q may be predicted from the induced stress σ and τ in the three dimensional direction (steps 106, 108).

Referring to FIG. 3a , the object 10 may be a wedge having two planar surfaces, a top surface 20 and an edge surface 30, each of which is subjected to normal load p on the planar surface and shear loads q along the planar surface respectively.

Referring to FIG. 3b , each of the normal load p and shear load q applied on the planar surfaces correspond to an exclusive equivalent normal p and shear q loads in half-space on the same planar surface. Such equivalent loads induce stresses σ, τ on the other plane surfaces of the wedge 10.

Turning now to the detailed numerical analysis of the present invention, each of the two surfaces 20, 30 of the wedge 10 may undergo discretization to form smaller, uniformly loaded sections. As shown in FIG. 4a , there is provided an effective grid meshing of the top surface 20 of the wedge 10. The surface 20 is divided into k rectangular grids and the ith rectangle 22 of which is shown in FIG. 4b . For instance, the ith rectangle may be dimensioned as 2α_(i) and 2β_(i) in the X₁ and Y₁ directions respectively. The coordinate of its centre is denoted as (x_(1i), y_(1i)) in the coordinate system X₁-Y₁-Z₁. The other surface 30 of the wedge 10 may also be discretized into a set of l rectangles each having a coordinate of (x_(2j), y_(2j)) in the coordinate system X₂-Y₂-Z₂.

In one example embodiment, the loads p, q applied on the two surfaces 20, 30 may be approached by piece-wise distributions. The values of the distributed loads on the two surfaces 20, 30 are represented by the values at the centres of rectangles. For instance, the distributed load of the ith rectangular section for the top surface 20 may be represented by (p₁)_(i)=p₁(x_(1i), y_(1i)), and (q_(y1))_(i)=q_(y1)(x_(1i), y_(1i)) and the distributed load of the jth rectangle for the edge surface 30 may be represented by (p₂)_(j)=p₂(x_(2j), y_(2j)), (q_(x2))_(j)=q_(x2)(x_(2j), y_(2j)) and (q_(y2))_(j)=q_(y2)(X_(2j), y_(2j)).

Preferably, the applied load applied on each of the plurality of surfaces 20 and 30 may be represented in the form of the accumulation of uniform loads applied on each uniformly loaded section. For instance, the normal and shear loads on the top surface 20 and the edge surface 30 may be defined as P₁, Q_(X1), Q_(Y1) and P₂, Q_(X2), Q_(Y2) in vector format, as expressed respectively in Equations (1) and (2)

$\begin{matrix} {{P_{1} = \begin{bmatrix} \left( p_{1} \right)_{1} \\ \left( p_{1} \right)_{2} \\ \vdots \\ \left( p_{1} \right)_{k} \end{bmatrix}};{Q_{X\; 1} = \begin{bmatrix} \left( q_{x\; 1} \right)_{1} \\ \left( q_{x\; 1} \right)_{2} \\ \vdots \\ \left( q_{x\; 1} \right)_{k} \end{bmatrix}};{Q_{Y\; 1} = \begin{bmatrix} \left( q_{y\; 1} \right)_{1} \\ \left( q_{y\; 1} \right)_{2} \\ \vdots \\ \left( q_{y\; 1} \right)_{k} \end{bmatrix}}} & (1) \\ {{P_{2} = \begin{bmatrix} \left( p_{2} \right)_{1} \\ \left( p_{2} \right)_{2} \\ \vdots \\ \left( p_{2} \right)_{l} \end{bmatrix}};{Q_{X\; 2} = \begin{bmatrix} \left( q_{x\; 2} \right)_{1} \\ \left( q_{x\; 2} \right)_{2} \\ \vdots \\ \left( q_{x\; 2} \right)_{l} \end{bmatrix}};{Q_{Y\; 2} = \begin{bmatrix} \left( q_{y\; 2} \right)_{1} \\ \left( q_{y\; 2} \right)_{2} \\ \vdots \\ \left( q_{y\; 2} \right)_{l} \end{bmatrix}}} & (2) \end{matrix}$

In one example embodiment, the distributed equivalent load of the ith rectangle for the top surface 20 may be represented by (p ₁)_(i)=p ₁(x_(1i), y_(1i)), (q _(x1))_(i)=q _(x1)(x_(1i), y_(1i)), (q _(y1))_(i)=q _(y1)(x_(1i), y_(1i)) and the distributed equivalent load of the jth rectangle for the edge surface 30 may be represented by (p ₂)_(i)=p ₂(x_(2j), y_(2j)), (q _(x2))_(j)=q _(x2)(x_(2j), y_(2j)), (q _(y2))_(j)=q _(y2)(x_(2j), y_(2j)).

In a similar manner, the equivalent load applied on each of the plurality of surfaces 20 and 30 may be represented in the form of the accumulation of uniform equivalent loads applied on each uniformly loaded section. For instance, the corresponding equivalent loads, P ₁, Q _(X1), and Q _(Y1), and P ₂, Q _(X2) and Q _(Y2) acting on the imaginary half-space formed by the extension of the top surface 20 of plane X₁-Y₁ and the edge surface 30 of plane X₂-Y₂ may also be expressed with piece-wise distribution in Equations (3) and (4).

$\begin{matrix} {{{\overset{\_}{P}}_{1} = \begin{bmatrix} \left( {\overset{\_}{p}}_{1} \right)_{1} \\ \left( {\overset{\_}{p}}_{1} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{p}}_{1} \right)_{k} \end{bmatrix}};{{\overset{\_}{Q}}_{X\; 1} = \begin{bmatrix} \left( {\overset{\_}{q}}_{x\; 1} \right)_{1} \\ \left( {\overset{\_}{q}}_{x\; 1} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{q}}_{x\; 1} \right)_{k} \end{bmatrix}};{{\overset{\_}{Q}}_{Y\; 1} = \begin{bmatrix} \left( {\overset{\_}{q}}_{y\; 1} \right)_{1} \\ \left( {\overset{\_}{q}}_{y\; 1} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{q}}_{y\; 1} \right)_{k} \end{bmatrix}}} & (3) \\ {{{\overset{\_}{P}}_{2} = \begin{bmatrix} \left( {\overset{\_}{p}}_{2} \right)_{1} \\ \left( {\overset{\_}{p}}_{2} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{p}}_{2} \right)_{l} \end{bmatrix}};{{\overset{\_}{Q}}_{X\; 2} = \begin{bmatrix} \left( {\overset{\_}{q}}_{x\; 2} \right)_{1} \\ \left( {\overset{\_}{q}}_{x\; 2} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{q}}_{x\; 2} \right)_{l} \end{bmatrix}};{{\overset{\_}{Q}}_{Y\; 2} = \begin{bmatrix} \left( {\overset{\_}{q}}_{y\; 2} \right)_{1} \\ \left( {\overset{\_}{q}}_{y\; 2} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{q}}_{y\; 2} \right)_{l} \end{bmatrix}}} & (4) \end{matrix}$

The inventors have further devised that the corresponding equivalent loads P ₁, Q _(X1), Q _(Y1) and P ₂, Q _(X2), Q _(Y2) are related to the original applied loads P₁, Q_(X1), Q_(Y1) and P₂, Q_(X2), Q_(Y2) onto the planar surface 20, 30 and a transformation matrix associated with the Poisson's ratio of the material, the wedge angle θ i.e. the orientation between the planar surfaces 20, 30. It may be appreciated by the persons skilled in the art that the equivalent load may be obtained from the applied load and the transformation matrix through mathematical calculations.

In one example embodiment, the equivalent load in vector format of each of the plurality of surfaces 20 and 30 may be represented by the transformation matrix containing two dimensional coefficients associated with the physical characteristic of the object 10 and the applied load. The transformation matrix may be derived from the accumulation of two dimensional coefficients of each uniformly loaded section used in the classical formulae for half-space stress and strain calculations. It would be appreciated by persons skilled in the art that the transformation matrix is dependent on the predetermined size and location of the uniformly loaded sections meshed from the plurality of surfaces 20 and 30.

In one example embodiment, the stress induced by equivalent loads on each of the plurality of planar surfaces 20 and 30 may also be represented in the form of the accumulation of stresses induced on the centre of each uniformly loaded section. For instance, the stresses at rectangle j on the plane X₂-Y₂ in the imaginary half-space induced by equivalent loads p ₁ and q ₁ on the plane X₁-Y₁ are denoted by their values at the centres of the rectangles, namely, (σ_(zz2))_(j)=σ_(zz2)(x_(2j), y_(2j)), (τ_(zx2))_(j)=σ_(zx2)(x_(2j), y_(2j), and (τ_(zy2))_(j)=τ_(zy2)(x_(2j), y_(2j)). These induced stresses of the plane X₂-Y₂ form three vectors, S_(Z), S_(X), and S_(Y) (Equation (5)).

$\begin{matrix} {{S_{Z} = \begin{bmatrix} \left( \sigma_{{zz}\; 2} \right)_{1} \\ \left( \sigma_{{zz}\; 2} \right)_{2} \\ \vdots \\ \left( \sigma_{{zz}\; 2} \right)_{l} \end{bmatrix}};{S_{X} = \begin{bmatrix} \left( \tau_{{zx}\; 2} \right)_{1} \\ \left( \tau_{{zx}\; 2} \right)_{2} \\ \vdots \\ \left( \tau_{{zx}\; 2} \right)_{l} \end{bmatrix}};{S_{Y} = \begin{bmatrix} \left( \tau_{{zy}\; 2} \right)_{1} \\ \left( \tau_{{zy}\; 2} \right)_{2} \\ \vdots \\ \left( \tau_{{zy}\; 2} \right)_{l} \end{bmatrix}}} & (5) \end{matrix}$

Similarly, the stresses at rectangle i on the plane X₁-Y₁ in the imaginary half-space induced by these equivalent loads on the plane X₂-Y₂ are denoted by their values at the centres of the rectangles, namely, (σ_(zz1))_(i)=σ_(zz1)(x_(1i), y_(1i)), (σ_(zx1))_(i)=τ_(zx1)(x_(1i), y_(1i)), (τ_(zy1))_(i)=τ_(zy1)(x_(1i), y_(1i)). These induced stresses of the plane X₁-Y₁ form three vectors, namely, T_(Z), T_(X), and T_(Y) (Equation (6)).

$\begin{matrix} {{T_{Z} = \begin{bmatrix} \left( \sigma_{{zz}\; 1} \right)_{1} \\ \left( \sigma_{{zz}\; 1} \right)_{2} \\ \vdots \\ \left( \sigma_{{zz}\; 1} \right)_{k} \end{bmatrix}};{T_{X} = \begin{bmatrix} \left( \tau_{{zx}\; 1} \right)_{1} \\ \left( \tau_{{zx}\; 1} \right)_{2} \\ \vdots \\ \left( \tau_{{zx}\; 1} \right)_{k} \end{bmatrix}};{T_{Y} = \begin{bmatrix} \left( \tau_{{zy}\; 1} \right)_{1} \\ \left( \tau_{{zy}\; 1} \right)_{2} \\ \vdots \\ \left( \tau_{{zy}\; 1} \right)_{k} \end{bmatrix}}} & (6) \end{matrix}$

The stresses σ_(zz2), τ_(zx2), and τ_(zy2) may be represented in the following expressions (Equations (7) to (9)) in terms of coefficients from m_(zzij) to m_(yyij) where the first subscript of m signifies the applied equivalent load on the ith rectangle on the plane X₁-Y₁. The second subscript represents the induced stress at the jth rectangle on the plane X₂-Y₂: σ_(zz2)(x _(2j,) y _(2j))=m _(zzij) p ₁(x _(1i) ,y _(1i))+m _(zxij) q _(x1)(x _(1i) ,y _(1i))+m _(zyij) q _(y1)(x _(1i) ,y _(1i))  (7) τ_(zx2)(x _(2j) ,y _(2j))=m _(xzij) p ₁(x _(1i) ,y _(1i))+m _(xxij) q _(x1)(x _(1i) ,y _(1i))+m _(xyij) q _(y1)(x _(1i) ,y _(1i))  (8) τ_(zy2)(x _(2j) ,y _(2j) =m _(yzij) p ₁(x _(1i) ,y _(1i))+m _(yxij) q _(x1)(x _(1i) ,y _(1i))+m _(yyij) q _(y1)(x _(1i) ,y _(1i))  (9)

Similarly, the stresses σ_(zz1), τ_(zx1), τ_(zy1) on the plane X₁-Y₁ induced by equivalent loads (p ₂)_(j), (q _(x2))_(j), (q _(y2))_(j) on the plane X₂-Y₂ are represented by the following expressions (Equations (10) to (12)): σ_(zz1)(x _(1i) ,y _(1i))n _(zzji) p ₂(x _(2j) ,y _(2j))+n _(zxji) q _(x2)(x _(2j) ,y _(2j))+n _(zyji) q _(y2)(x _(2j) ,y _(2j))  (10) τ_(zx1)(x _(1i) ,y _(1i))=n _(xzji) p ₂(x _(2j) ,y _(2j))+n _(xxji) q _(x2)(x _(2j) ,y _(2j))+n _(xyji) q _(y2)(x _(2j) ,y _(2j))  (11) τ_(zy1)(x _(1i) ,y _(1i))=n _(yzji) p ₂(x _(2j) ,y _(2j))+n _(yxji) q _(x2)(x _(2i) ,y _(2j))+n _(yyji) q _(y2)(x _(2j) ,y _(2j))  (12)

The relationship between these equivalent loads and stresses of half-space model may be determined. Accordingly, the three induced stress vectors S_(Z), S_(X) S_(Y) on the plane X₂-Y₂ due to the equivalent normal and shear load vectors P ₁, Q _(X1), Q _(Y1) on the plane X₁-Y₁ can be determined based on the half-space model by the following expressions: S _(Z) =M _(ZZ) ·P ₁ +M _(ZX) ·Q _(X1) +M _(ZY) ·Q _(Y1)  (13) S _(X) =M _(XZ) ·P ₁ +M _(XX) ·Q _(X1) +M _(XY) ·Q _(Y1)  (14) S _(Y) =M _(YZ) ·P ₁ +M _(YX) ·Q _(X1) +M _(YY) ·Q _(Y1)  (15)

Similarly, the equivalent normal and shear load vectors, P ₂, Q _(X2), Q _(Y2) on the surface X₂-Y₂ of half-space induces stress vectors, T_(Z), T_(X), and T_(Y), on the plane X₁-Y₁ can be determined based on the half-space model by: T _(Z) =N _(ZZ) ·P ₂ +N _(ZX) ·Q _(X2) +N _(ZY) ·Q _(Y2)  (16) T _(X) =N _(XZ) ·P ₂ +N _(XX) ·Q _(X2) +N _(XY) ·Q _(Y2)  (17) T _(Y) =N _(YZ) ·P ₂ +N _(YX) ·Q _(X2) +N _(YY) ·Q _(Y2)  (18)

The two-dimensional array of the coefficients forms the following nine reflecting matrices M_(ZZ), M_(ZX), M_(ZY), M_(XZ), M_(XX), M_(XY), M_(YZ), M_(YX) and M_(Y) (in Equations (19) to (27)) and N_(ZZ), N_(ZX), N_(ZY), N_(XZ), N_(XX), N_(XY), N_(YZ), N_(YX) and N_(YY) (in Equations (28) to (36)).

$\begin{matrix} {M_{ZZ} = \begin{bmatrix} m_{{zz}\; 11} & m_{{zz}\; 21} & \ldots & m_{{zzk}\; 1} \\ m_{{zz}\; 12} & m_{{zz}\; 22} & \ldots & m_{{zzk}\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ m_{{zz}\; 1l} & m_{{zz}\; 2l} & \ldots & m_{zzkl} \end{bmatrix}} & (19) \\ {M_{ZX} = \begin{bmatrix} m_{{zx}\; 11} & m_{{zx}\; 21} & \ldots & m_{{zxk}\; 1} \\ m_{{zx}\; 12} & m_{{zx}\; 22} & \ldots & m_{{zxk}\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ m_{{zx}\; 1l} & m_{{zx}\; 2l} & \ldots & m_{zxkl} \end{bmatrix}} & (20) \\ \vdots & \vdots \\ {M_{YY} = \begin{bmatrix} m_{{yy}\; 11} & m_{{yy}\; 21} & \ldots & m_{{yyk}\; 1} \\ m_{{yy}\; 12} & m_{{yy}\; 22} & \ldots & m_{{yyk}\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ m_{{yy}\; 1l} & m_{{yy}\; 2l} & \ldots & m_{yykl} \end{bmatrix}} & (27) \\ {N_{ZZ} = \begin{bmatrix} n_{{zz}\; 11} & n_{{zz}\; 21} & \ldots & n_{{zzl}\; 1} \\ n_{{zz}\; 12} & n_{{zz}\; 22} & \ldots & n_{{zzl}\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ n_{{zz}\; 1k} & n_{{zz}\; 2k} & \ldots & n_{zzlk} \end{bmatrix}} & (28) \\ {N_{ZX} = \begin{bmatrix} n_{{zx}\; 11} & n_{{zx}\; 21} & \ldots & n_{{zxl}\; 1} \\ n_{{zx}\; 12} & n_{{zx}\; 22} & \ldots & n_{{zxl}\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ n_{{zx}\; 1k} & n_{{zx}\; 2k} & \ldots & n_{zxlk} \end{bmatrix}} & (29) \\ \vdots & \vdots \\ {N_{YY} = \begin{bmatrix} n_{{yy}\; 11} & n_{{yy}\; 21} & \ldots & n_{{yyl}\; 1} \\ n_{{yy}\; 12} & n_{{yy}\; 22} & \ldots & n_{{yyl}\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ n_{{yy}\; 1k} & n_{{yy}\; 2k} & \ldots & n_{yylk} \end{bmatrix}} & (36) \end{matrix}$

By knowing the physical properties of the wedge 10 and the normal and shear loads P₁, Q_(X1), Q_(Y1) and P₂, Q_(X2), Q_(Y2) applied on the top and edge surfaces 20, 30, the equivalent loads P ₁, Q _(X1), and Q _(Y1), and P ₂, Q _(X2) and Q _(Y2) on the two wedge surfaces 20 and 30 may be determined readily from the transformation matrix by the following expression:

$\begin{matrix} {\begin{bmatrix} {\overset{\_}{P\;}}_{1} \\ {\overset{\_}{Q}\;}_{X\; 1} \\ {\overset{\_}{Q}\;}_{Y\; 1} \\ {{\overset{\_}{P\;}}_{2}\;} \\ {\overset{\_}{Q}\;}_{X\; 2} \\ {\overset{\_}{Q}\;}_{Y\; 2} \end{bmatrix} = {\begin{bmatrix} I & 0 & 0 & {- N_{ZZ}} & {- N_{ZX}} & {- N_{ZY}} \\ {0\;} & I & 0 & {- N_{XZ}} & {- N_{XX}} & {- N_{XY}} \\ 0 & 0 & I & {- N_{YZ}} & {- N_{YX}} & {- N_{YY}} \\ {- M_{ZZ}} & {- M_{ZX}} & {- M_{ZY}} & I & 0 & 0 \\ {- M_{XZ}} & {- M_{XX}} & {- M_{XY}} & 0 & I & 0 \\ {- M_{YZ}} & {- M_{YX}} & {- M_{YY}} & 0 & 0 & I \end{bmatrix}^{- 1}\begin{bmatrix} P_{1} \\ Q_{X\; 1} \\ Q_{Y\; 1} \\ P_{2} \\ Q_{X\; 2} \\ Q_{Y\; 2} \end{bmatrix}}} & (37) \end{matrix}$ The transformation matrix is expressed with:

$\begin{matrix} \begin{bmatrix} I & 0 & 0 & {- N_{ZZ}} & {- N_{ZX}} & {- N_{ZY}} \\ {0\;} & I & 0 & {- N_{XZ}} & {- N_{XX}} & {- N_{XY}} \\ 0 & 0 & I & {- N_{YZ}} & {- N_{YX}} & {- N_{YY}} \\ {- M_{ZZ}} & {- M_{ZX}} & {- M_{ZY}} & I & 0 & 0 \\ {- M_{XZ}} & {- M_{XX}} & {- M_{XY}} & 0 & I & 0 \\ {- M_{YZ}} & {- M_{YX}} & {- M_{YY}} & 0 & 0 & I \end{bmatrix}^{- 1} & (38) \end{matrix}$

By determining the equivalent load P ₁, Q _(X1), and Q _(Y1), and P ₂, Q _(X2) and Q _(Y2) on the two edge surfaces 20 and 30 with the transformation matrix above (Equation (38)), the induced stress vectors may be obtained from the known equivalent loads using the classical half-space formulae accordingly.

In one example embodiment, the analytical method 100 in the present invention is adopted for determining the stress and strain of a wedge 10 with a Poisson's ratio of 0.3, a load applied on the top surface 20 of the wedge 10, and a wedge angle θ between the top surface 20 and the edge surface 30 of the wedge 10, as shown in FIG. 5. The top surface 20 is divided into 8500 rectangular grids with 85 units in the Y₁-direction and 100 units in the X₁-direction.

For instance, the load applied on the top surface 20 may be a half-semispherical normal load distribution 40, as schematically shown in FIGS. 5 and 6 and represented in the following Equation (39). The load is applied to the top surface 20 (Z1=0) of the elastic wedge 10 with the maximum, p₀, at the origin (0, 0, 0).

$\begin{matrix} {{{{p\left( {x_{1},y_{1}} \right)} = {p_{0}\frac{\sqrt{a^{2} - x_{1}^{2} - y_{1}^{2}}}{a}}},{\sqrt{x_{1}^{2} + y_{1}^{2}} \leq a}}{{{p\left( {x_{1},y_{1}} \right)} = 0},{\sqrt{x_{1}^{2} + y_{1}^{2}} > a}}} & (39) \end{matrix}$ where a is the radius of the half-semispherical load distribution.

In one specific embodiment, the two planar surfaces 20 and 30 forms a wedge angle θ of 170°. The six equivalent loads on the top planar surface 20 (as shown in FIGS. 7a to 7c ) and edge planar surface 30 (as shown in FIGS. 7d to 7f ) are determined from the original load and the transformation matrix (step 104). As shown in FIG. 7a , the equivalent load 50 which is in the direction normal to the top surface 20 of the wedge 10 deviates very little from the original applied load 40 as shown in FIG. 6. The other equivalent loads as shown in FIGS. 7b to 7f are of much smaller magnitude. The resultant stress is thus dominated by the equivalent load 50.

The internal stress field of the wedge 10 is determined from the summation of stresses calculated with the individual equivalent load (step 106), and the normal and shear stresses on y=0 plane are presented in FIGS. 8a to 8c . Finally, the mechanical failure of the wedge 10 is predicted in FIG. 8d (step 108).

In another specific embodiment, the two planar surfaces 20 and 30 forms a wedge angle θ of 90°. The six equivalent loads on the two planar surfaces 20 and 30 are determined from the original load and the transformation matrix (step 104). The normal and shear stresses are presented in FIGS. 9a to 9c (step 106), and thus the von Mises stress is determined, as shown in FIG. 9d . Similarly, the strain field can be calculated (step 108). Thus, mechanical failure of the wedge 10 is predicted.

In yet another specific embodiment, the two planar surfaces 20 and 30 forms a wedge angle θ of 60°. The six equivalent loads on the two planar surfaces 20 and 30 are determined from the original load and the transformation matrix (step 104). The normal and shear stresses on y=0 plane are presented in FIGS. 10a to 10c (step 106), and thus the von Mises stress is determined, as shown in FIG. 10d . Similarly, the strain field can be calculated (step 108). Thus, mechanical failure of the wedge 10 is predicted.

In one further example embodiment, the stress behaviour of the object 10 with a wedge angle θ of 135° and 70° are also determined by using the present invention and each presented in FIGS. 11a to 11c and FIGS. 12a to 12c respectively.

As shown in FIGS. 8a to 12c , it shows the difference in the normal stresses in the y=0 plane with varying wedge angles θ. The maximum p_(o) of the half-semispherical distributed-load, as expressed in Equation (35), is at the origin (0, 0, 0).

Referring to the example of σ_(xx) in FIGS. 8a, 9a, 10a, 11a and 12a , the maximum tensile stress located on the top surface 20 next to the right boundary of the applied load increases with decreasing wedge angle θ. The maximum σ_(xx) for the wedge with 170° is 0.089 p₀ on the top surface, but the maximum σ_(xx) reaches 0.41 p₀ for quarter-space. The compressive stress on the top surface 20 decreases with the wedge angle θ. Specially, the compressive stress on the top surface 20 totally changes into tensile stress when the wedge angle θ is <90°. As it would be appreciated by person skilled in the art, this finding would be reasonable given the lack of enough support under the top surface 20 for such wedge structure 10 with angle θ<90°.

Turning now to the normal stress σ_(xx) inside the contact body 10, the trend can be separated into two parts. For wedges with angle θ>90°, the area of compressive stress shrinks to the apex of the wedge 10 with decreasing wedge angle θ. However, the area of compressive stress increases under the load area with the wedge angle θ decreasing from 90°. Furthermore, the compressive stress along the free surface decreases slightly with depth for the wedge with 60°.

The stresses σ_(yy) are totally compressive when the wedge angle θ is >90°. However, tensile stress appears next to the free surface when the wedge 10 is small enough, such as the wedge θ with 60°. Additionally, the effect of wedge angle θ on σ_(yy) in the corner is apparent. The gradient of σ_(yy) at the corner drops gradually with the decrease in wedge angle θ from 170° to 90°. However, this trend is reversed once the wedge angle θ reaches below 90°. A similar finding was observed for normal stress σ_(zz).

In one specific example embodiment, the analytical method 100 in the present invention is adopted for determining the stress and strain of an object 10 with two identical loads being applied onto the top 20 and edge surfaces 30 with a wedge angle θ there between simultaneously.

For instance, a pair of normal load distributions, represented in the following Equation (40), is applied to the top surface 20 and the edge surface 30 respectively:

$\begin{matrix} {{{{p\left( {x_{1},y_{1}} \right)} = {p_{0}\frac{\sqrt{a^{2} - \left( {x_{1} - a} \right)^{2} - y_{1}^{2}}}{a}}},{\sqrt{\left( {x_{1} - a} \right)^{2} + y_{1}^{2}} \leq a}}{{{p\left( {x_{1},y_{1}} \right)} = 0},{\sqrt{\left( {x_{1} - a} \right)^{2} + y_{1}^{2}} > a}}} & (40) \end{matrix}$

The normal and shear stresses on y=0 plane are presented in FIGS. 13a to 13b respectively (step 106). In detail, the maximum modulus of σ_(xx)/p₀ is 1.0 and just locates at (0, 0, 1.0), which is the centre of the side semispherical surface load.

As the current wedge problem is symmetrical about the plane x1-z1, it would be appreciated by person skilled in the art that the result of σ_(zz) would also be depicted by FIG. 13(a). The distribution of σ_(yy) in y=0 plane is symmetrical about x1-z1 (as shown in FIG. 13b ). Accordingly, the maximum stress of the object 10 is determined to be present at (0.5, 0, 0.5), as shown in FIG. 13 c.

Advantageously, the method in the present invention requires only the two side surfaces 20 and 30 of the wedge structure 10 is meshed into rectangular grids. Therefore, it is much simpler and more suitable to study the wedge problems with varying surface loads.

Furthermore, the transformation matrix provide in the present invention is only related to the wedge angle θ, mesh structure and the elastic properties of object 10 but not the applied load p. The transformation matrix can be determined and saved in advance for subsequent calculations with different loads p, q. Therefore, the explicit matrix method is suitable to be embedded into those studies requiring iterative calculations of wedge deformation or stress, such as the elastohydrodynamic lubrication analysis of roller bearings.

It will be appreciated by persons skilled in the art that numerous variations and/or modifications may be made to the invention as shown in the specific embodiments without departing from the spirit or scope of the invention as broadly described. The present embodiments are, therefore, to be considered in all respects as illustrative and not restrictive.

Any reference to prior art contained herein is not to be taken as an admission that the information is common general knowledge, unless otherwise indicated. 

The invention claimed is:
 1. An analytical method of determining the stress and strain of an object with a plurality of planar surfaces under one or more applied loads, comprising the steps of: A. applying three dimensional loads on each of the plurality of planar surfaces simultaneously; B. determining the three dimensional equivalent loads applied on each of the plurality of planar surfaces from the applied three dimensional loads with the stress boundary conditions; C. determining the stress and strain fields induced on each of the plurality of planar surfaces by the equivalent load applied on the other plurality of planar surfaces; and D. determining the stress and strain associated with the mechanical failure of the object; wherein step B includes determining a transformation matrix associated with the physical characteristics of the object and in turn the equivalent load from the applied load; wherein each of the plurality of planar surfaces is divided into a plurality of uniformly loaded sections each having a centre respectively; wherein each centre of the uniformly loaded sections on the plurality of planar surfaces X₁-Y₁ and X₂-Y₂ is represented by x_(1i), y_(1i) and x_(2i), y_(2i) respectively; and wherein the parameter of the applied load P₁, Q_(x1), Q_(y1) applied on each of the plurality of planar surfaces Xi-Y₁ and the parameter of the applied load p₂, Q_(x2), Q_(y2), applied on each of the plurality of planar surfaces X₂-Y₂ are each represented in the form of the accumulation of uniform loads p, q applied on each uniformly loaded section by the following expression respectively: ${P_{1} = \begin{bmatrix} \left( p_{1} \right)_{1} \\ \left( p_{1} \right)_{2} \\ \vdots \\ \left( p_{1} \right)_{k} \end{bmatrix}};{Q_{X\; 1} = \begin{bmatrix} \left( q_{x\; 1} \right)_{1} \\ \left( q_{x\; 1} \right)_{2} \\ \vdots \\ \left( q_{x\; 1} \right)_{k} \end{bmatrix}};{Q_{Y\; 2} = \begin{bmatrix} \left( q_{y\; 1} \right)_{1} \\ \left( q_{y\; 1} \right)_{2} \\ \vdots \\ \left( q_{y\; 1} \right)_{k} \end{bmatrix}}$ ${P_{2} = \begin{bmatrix} \left( p_{2} \right)_{1} \\ \left( p_{2} \right)_{2} \\ \vdots \\ \left( p_{2} \right)_{l} \end{bmatrix}};{Q_{X\; 2} = \begin{bmatrix} \left( q_{x\; 2} \right)_{1} \\ \left( q_{x\; 2} \right)_{2} \\ \vdots \\ \left( q_{x\; 2} \right)_{l} \end{bmatrix}};{Q_{Y\; 2} = \begin{bmatrix} \left( q_{y\; 2} \right)_{1} \\ \left( q_{y\; 2} \right)_{2} \\ \vdots \\ \left( q_{y\; 2} \right)_{l} \end{bmatrix}}$ where the value of each uniform load p, q on the uniformly loaded section is represented by the value at the center (x_(1i), y_(1i)), (x₂, y₂) of the section respectively: (p₁)_(i)=p₁(x_(1i), y_(1i)), (q_(x1))_(i)=q_(x1)(x_(1i), y_(1i)), (q_(y1))=q_(y1)(x_(1i), y_(1i)); and (p₂)_(i)=p₂(x_(2i), y_(2i)), (q_(x2))_(i)=q_(x2)(x_(2i), y_(2i))_(i)=q_(y2)(x_(2j), y_(2j)).
 2. The analytical method according to claim 1, wherein the applied load includes the accumulation of uniform loads applied on each uniformly loaded section.
 3. The analytical method according to claim 1, wherein the equivalent load includes the accumulation of uniform equivalent loads applied on each uniformly loaded section.
 4. The analytical method according to claim 1, wherein the induced stress includes the accumulation of stresses induced on the centre of each uniformly loaded section.
 5. The analytical method according to claim 1, wherein the transformation matrix is associated with the predetermined size of the uniformly loaded section.
 6. An analytical method of determining the stress and strain of an object with a plurality of planar surfaces under one or more applied loads, comprising the steps of: A. applying three dimensional loads on each of the plurality of planar surfaces simultaneously; B. determining the three dimensional equivalent loads applied on each of the plurality of planar surfaces from the applied three dimensional loads with the stress boundary conditions; C. determining the stress and strain fields induced on each of the plurality of planar surfaces by the equivalent load applied on the other plurality of planar surfaces; and D. determining the stress and strain associated with the mechanical failure of the object; wherein step B includes determining a transformation matrix associated with the physical characteristics of the object and in turn the equivalent load from the applied load; wherein each of the plurality of planar surfaces is divided into a plurality of uniformly loaded sections each having a centre respectively; wherein each centre of the uniformly loaded sections on the plurality of planar surfaces X₁-Y₁ and X₂-Y₂ is represented by x_(1i), y_(1i) and x_(2i), y_(2i) respectively; and wherein the parameter of the equivalent load P ₁, Q _(X1), Q _(Y1) applied on the plurality of planar surfaces X₁-Y₁ and the parameter of the equivalent load P ₂, Q _(X2), Q _(Y2) applied on the plurality of planar surfaces X₂-Y₂ are each represented in the form of the accumulation of equivalent loads P, Q applied on each of uniformly loaded section by the following expression respectively: ${{\overset{\_}{P}}_{1} = \begin{bmatrix} \left( {\overset{\_}{p}}_{1} \right)_{1} \\ \left( {\overset{\_}{p}}_{1} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{p}}_{1} \right)_{k} \end{bmatrix}};{{\overset{\_}{Q}}_{X\; 1} = \begin{bmatrix} \left( {\overset{\_}{q}}_{x\; 1} \right)_{1} \\ \left( {\overset{\_}{q}}_{x\; 1} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{q}}_{x\; 1} \right)_{k} \end{bmatrix}};{{\overset{\_}{Q}}_{Y\; 1} = \begin{bmatrix} \left( {\overset{\_}{q}}_{y\; 1} \right)_{1} \\ \left( {\overset{\_}{q}}_{y\; 1} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{q}}_{y\; 1} \right)_{k} \end{bmatrix}}$ ${{\overset{\_}{P}}_{2} = \begin{bmatrix} \left( {\overset{\_}{p}}_{2} \right)_{1} \\ \left( {\overset{\_}{p}}_{2} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{p}}_{2} \right)_{l} \end{bmatrix}};{{\overset{\_}{Q}}_{X\; 2} = \begin{bmatrix} \left( {\overset{\_}{q}}_{x\; 2} \right)_{1} \\ \left( {\overset{\_}{q}}_{x\; 2} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{q}}_{x\; 2} \right)_{l} \end{bmatrix}};{{\overset{\_}{Q}}_{Y\; 2} = \begin{bmatrix} \left( {\overset{\_}{q}}_{y\; 2} \right)_{1} \\ \left( {\overset{\_}{q}}_{y\; 2} \right)_{2} \\ \vdots \\ \left( {\overset{\_}{q}}_{y\; 2} \right)_{l} \end{bmatrix}}$ where the value of each uniform equivalent load p, q on the uniformly loaded section is represented by the value at the center (x_(1i)y_(1i)), (x₂y₂) of the section respectively: (p ₁)_(i)=p ₁(x_(1i), y_(1i)), (q _(x1))_(i)=q _(x1)(x_(1i), y_(1i)), (q _(y1))_(i)=q _(y1)(x_(1i), y_(1i)); and (p ₂)_(i)=p ₂(x_(2i), y_(2i)), (q _(x2))_(i)=q _(x2)(x_(2i), y_(2i)), (q _(y2))_(i)=q _(y2)(x_(2j), y_(2j)).
 7. The analytical method according to claim 6, wherein the applied load includes the accumulation of uniform loads applied on each uniformly loaded section.
 8. The analytical method according to claim 6, wherein the equivalent load includes the accumulation of uniform equivalent loads applied on each uniformly loaded section.
 9. The analytical method according to claim 6, wherein the induced stress includes the accumulation of stresses induced on the centre of each uniformly loaded section.
 10. The analytical method according to claim 6, wherein the transformation matrix is associated with the predetermined size of the uniformly loaded section.
 11. The analytical method according to claim 6, wherein equivalent load P ₁, Q _(X1), Q _(Y1) and equivalent load P ₂, Q _(X2), Q _(Y2) are determined from the following expression: $\begin{bmatrix} {\overset{\_}{P}\;}_{1} \\ {\overset{\_}{Q}\;}_{X\; 1} \\ {\overset{\_}{Q}\;}_{Y\; 1} \\ {\overset{\_}{P}\;}_{2} \\ {\overset{\_}{Q}\;}_{X\; 2} \\ {\overset{\_}{Q}\;}_{Y\; 2} \end{bmatrix} = {\begin{bmatrix} I & 0 & 0 & {- N_{ZZ}} & {- N_{ZX}} & {- N_{ZY}} \\ {0\;} & I & 0 & {- N_{XZ}} & {- N_{XX}} & {- N_{XY}} \\ 0 & 0 & I & {- N_{YZ}} & {- N_{YX}} & {- N_{YY}} \\ {- M_{ZZ}} & {- M_{ZX}} & {- M_{ZY}} & I & 0 & 0 \\ {- M_{XZ}} & {- M_{XX}} & {- M_{XY}} & 0 & I & 0 \\ {- M_{YZ}} & {- M_{YX}} & {- M_{YY}} & 0 & 0 & I \end{bmatrix}^{- 1}\begin{bmatrix} P_{1} \\ Q_{X\; 1} \\ Q_{Y\; 1} \\ P_{2} \\ Q_{X\; 2} \\ Q_{Y\; 2} \end{bmatrix}}$ where the inverse matrix is the transformation matrix; M_(ZZ), M_(ZX), M_(ZY), M_(XZ), M_(XX), M_(XY), M_(YZ), M_(YX)M_(YY) are the two dimensional matrices of coefficients for stress calculation using classical half-space formulae based on equivalent loads on planar surface X₁-Y₁; N_(ZZ), N_(ZX)N_(ZY), N_(XZ), N_(XX)N_(XY), N_(YZ), N_(YX), N_(YY) are the two dimensional matrices of coefficients for stress calculation using classical half-space formulae based on equivalent loads on planar surface X₂-Y₂; P₁, Q_(X1), Q_(Y1) are the loads applied on planar surface X₁-Y₁; and P_(Z), Q_(XZ), Q_(YZ) are the loads applied on planar surface X₂-Y₂.
 12. An analytical method of determining the stress and strain of an object with a plurality of planar surfaces under one or more applied loads, comprising the steps of: A. applying three dimensional loads on each of the plurality of planar surfaces simultaneously; B. determining the three dimensional equivalent loads applied on each of the plurality of planar surfaces from the applied three dimensional loads with the stress boundary conditions; C. determining the stress and strain fields induced on each of the plurality of planar surfaces by the equivalent load applied on the other plurality of planar surfaces; and D. determining the stress and strain associated with the mechanical failure of the object; wherein step B includes determining a transformation matrix associated with the physical characteristics of the object and in turn the equivalent load from the applied load; and wherein the transformation matrix is represented by a plurality of two dimensional coefficients M_(ZZ), M_(ZX), M_(ZY), M_(XZ), M_(XX), M_(YZ), M_(YX), M_(YY) and N_(ZZ), N_(ZX), N_(ZY), N_(XZ), N_(xx)N_(XY), N_(YZ)N_(YX), N_(YY) respectively in the form of $\begin{bmatrix} I & 0 & 0 & {- N_{ZZ}} & {- N_{ZX}} & {- N_{ZY}} \\ {0\;} & I & 0 & {- N_{XZ}} & {- N_{XX}} & {- N_{XY}} \\ 0 & 0 & I & {- N_{YZ}} & {- N_{YX}} & {- N_{YY}} \\ {- M_{ZZ}} & {- M_{ZX}} & {- M_{ZY}} & I & 0 & 0 \\ {- M_{XZ}} & {- M_{XX}} & {- M_{XY}} & 0 & I & 0 \\ {- M_{YZ}} & {- M_{YX}} & {- M_{YY}} & 0 & 0 & I \end{bmatrix}^{- 1},{and}$ wherein M_(ZZ), M_(ZX), M_(ZY), M_(XZ), M_(XX), M_(XY), M_(YZ), M_(YX), M_(YY) and N_(ZZ), N_(ZX), N_(ZY), N_(XZ), N_(XX), N_(XY), N_(YZ), N_(YX), N_(YY) are each represented in the form of the accumulation of coefficients of each uniformly loaded section by the following expression respectively: $M_{ZZ} = \begin{bmatrix} m_{{zz}\; 11} & m_{{zz}\; 21} & \ldots & m_{{zzk}\; 1} \\ m_{{zz}\; 12} & m_{{zz}\; 22} & \ldots & m_{{zzk}\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ m_{{zz}\; 1l} & m_{{zz}\; 2l} & \ldots & m_{zzkl} \end{bmatrix}$ $\;{M_{ZX} = \begin{bmatrix} m_{{zx}\; 11} & m_{{zx}\; 21} & \ldots & m_{{zxk}\; 1} \\ m_{{zx}\; 12} & m_{{zx}\; 22} & \ldots & m_{{zxk}\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ m_{{zx}\; 1l} & m_{{zx}\; 2l} & \ldots & m_{zxkl} \end{bmatrix}}$ ⋮ $M_{YY} = \begin{bmatrix} m_{{yy}\; 11} & m_{{yy}\; 21} & \ldots & m_{{yyk}\; 1} \\ m_{{yy}\; 12} & m_{{yy}\; 22} & \ldots & m_{{yyk}\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ m_{{yy}\; 1l} & m_{{yy}\; 2l} & \ldots & m_{yykl} \end{bmatrix}$ $N_{ZZ} = \begin{bmatrix} n_{{zz}\; 11} & n_{{zz}\; 21} & \ldots & n_{{zzl}\; 1} \\ n_{{zz}\; 12} & n_{{zz}\; 22} & \ldots & n_{{zzl}\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ n_{{zz}\; 1k} & n_{{zz}\; 2k} & \ldots & n_{zzlk} \end{bmatrix}$ $N_{ZX} = \begin{bmatrix} n_{{zx}\; 11} & n_{{zx}\; 21} & \ldots & n_{{zxl}\; 1} \\ n_{{zx}\; 12} & n_{{zx}\; 22} & \ldots & n_{{zxl}\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ n_{{zx}\; 1k} & n_{{zx}\; 2k} & \ldots & n_{zxlk} \end{bmatrix}$ ⋮ $N_{YY} = {\begin{bmatrix} n_{{yy}\; 11} & n_{{yy}\; 21} & \ldots & n_{{yyl}\; 1} \\ n_{{yy}\; 12} & n_{{yy}\; 22} & \ldots & n_{{yyl}\; 2} \\ \vdots & \vdots & \ddots & \vdots \\ n_{{yy}\; 1k} & n_{{yy}\; 2k} & \ldots & n_{yylk} \end{bmatrix}.}$
 13. The analytical method according to claim 12, wherein the physical characteristic includes the elastic property of a said object and the relative orientations between the planar surfaces.
 14. The analytical method according to claim 12, wherein the transformation matrix includes the accumulation of coefficients used in the classical formulae for half-space stress and strain calculations.
 15. The analytical method according to claim 12, wherein the parameter of the induced stress S_(z), S_(X), S_(Y) induced by the equivalent load P _(Z), Q _(XZ), Q _(YZ) on the plurality of planar surfaces X₂-Y₂ and the parameter of the induced stress T_(Z), T_(X), T_(Y) induced by the equivalent load P ₁, Q _(X1), Q _(Y1) on the plurality of planar surfaces X₁-Y₁ are each represented in the form of the accumulation of stresses σ, τ induced on the centre of each uniformly loaded section by the following expression respectively: ${S_{Z} = \begin{bmatrix} \left( \sigma_{{zz}\; 2} \right)_{1} \\ \left( \sigma_{{zz}\; 2} \right)_{2} \\ \vdots \\ \left( \sigma_{{zz}\; 2} \right)_{l} \end{bmatrix}};{S_{X} = \begin{bmatrix} \left( \tau_{{zx}\; 2} \right)_{1} \\ \left( \tau_{{zx}\; 2} \right)_{2} \\ \vdots \\ \left( \tau_{{zx}\; 2} \right)_{l} \end{bmatrix}};{S_{Y} = \begin{bmatrix} \left( \tau_{{zy}\; 2} \right)_{1} \\ \left( \tau_{{zy}\; 2} \right)_{2} \\ \vdots \\ \left( \tau_{{zy}\; 2} \right)_{l} \end{bmatrix}}$ ${T_{Z} = \begin{bmatrix} \left( \sigma_{{zz}\; 1} \right)_{1} \\ \left( \sigma_{{zz}\; 1} \right)_{2} \\ \vdots \\ \left( \sigma_{{zz}\; 1} \right)_{k} \end{bmatrix}};{T_{X} = \begin{bmatrix} \left( \tau_{{zx}\; 1} \right)_{1} \\ \left( \tau_{{zx}\; 1} \right)_{2} \\ \vdots \\ \left( \tau_{{zx}\; 1} \right)_{k} \end{bmatrix}};{T_{Y} = \begin{bmatrix} \left( \tau_{{zy}\; 1} \right)_{1} \\ \left( \tau_{{zy}\; 1} \right)_{2} \\ \vdots \\ \left( \tau_{{zy}\; 1} \right)_{k} \end{bmatrix}}$ where the value of each stresses σ, τ on the centre of each uniformly loaded section is represented by the value at the center (x_(1i), y_(1i)), (x₂, y₂ of the section respectively: (σ_(zz2))_(i)=σ_(zz2)(x_(2i), y_(2i)), (τ_(zx2))_(i)=τ_(zx2)(x_(2i), y_(2i)), (τ_(zy2))_(i)=τ_(zy2)(x_(2j), y_(2j)); and (σ_(zz1))_(i)=σ_(zz1)(x_(1i), y_(1i)), (τ_(zx1))_(i)=τ_(zx1)(x_(1i), y_(1i)), (τ_(zy1))_(i)=τ_(zy1)(x_(1i), y_(1i)).
 16. The analytical method according to claim 15, wherein the parameter of the induced stress S_(Z), S_(X), S_(Y) induced on the plurality of planar surfaces X₂-Y₂ and the parameter of the induced stress T_(Z), T_(Z), T_(Y) induced on the plurality of planar surfaces X₁-Y₁ are each determined from the following expressions respectively: S _(Z) =M _(ZZ) ·P ₁ +M _(ZX) ·Q _(X1) +M _(ZY) ·Q _(Y1) S _(X) =M _(XZ) ·P ₁ +M _(XX) ·Q _(X1) +M _(XY) ·Q _(Y1) S _(Y) =M _(YZ) ·P ₁ +M _(YX) ·Q _(X1) +M _(YY) ·Q _(Y1) T _(Z) =N _(ZZ) ·P ₂ +N _(ZX) ·Q _(X2) +N _(ZY) ·Q _(Y2) T _(X) =N _(XZ) ·P ₂ +N _(XX) ·Q _(X2) +N _(XY) ·Q _(Y2) T _(Y) =N _(YZ) ·P ₂ +N _(YX) ·Q _(X2) +N _(YY) ·Q _(Y2) where M_(ZZ), M_(ZX), M_(XY), M_(XZ), M_(XX), M_(YY), M_(YZ), M_(YX), M_(YY) and N_(ZZ), N_(ZX), N_(ZY), N_(XZ), N_(XX), N_(XY), N_(YZ), N_(YX), N_(YY) are two dimensional coefficients. 